A Class of Special Matrices and Quantum Entanglement 



Shao-Ming Fei^ and Xianqing Li-Jost* 

^Department of Mathematics, Capital Normal University, Beijing 100037, P.R. China 
^Institut fur Angewandte Mathematik, Universitat Bonn, 53115 Bonn, Germany 
E-mail : f ei @ wiener . iam . uni-b onn . de 

*Max-Planck-Institute for Mathematics in the Sciences, 04103 Leipzig, Germany 
E-mail: Xianqing.Li-Jost@mis.mpg.de 



Abstract 

We present a kind of construction for a class of special matrices with 
at most two different eigenvalues, in terms of some interesting multi- 
plicators which are very useful in calculating eigenvalue polynomials of 
these matrices. This class of matrices defines a special kind of quan- 
tum states — (^-computable states. The entanglement of formation for 
a large class of quantum mixed states is explicitly presented. 
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1 Introduction 

Quantum entangled states are playing an important role in quantum communication, infor- 
mation processing and quantum computing [1], especially in the investigation of quantum 
teleportation [2, 3], dense coding [5], decoherence in quantum computers and the evaluation 
of quantum cryptographic schemes [6]. To quantify entanglement, a number of entanglement 
measures such as the entanglement of formation and distillation [7, 8, 9], negativity [10, 11], 
relative entropy [9, 12] have been proposed for bipartite states [6, 8] [11-13]. Most of these 
measures of entanglement involve extremizations which are difficult to handle analytically. 
For instance the entanglement of formation [7] is intended to quantify the amount of quan- 
tum communication required to create a given state. The entanglement of formation for a 
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pair of qubits can be expressed as a monotonically increasing function of the "concurrence" , 
which can be taken as a measure of entanglement in its own right [14]. From the expression 
of this concurrence, the entanglement of formation for mixed states of a pair of qubits is 
calculated [14]. Although entanglement of formation is defined for arbitrary dimension, so 
far no explicit analytic formulae for entanglement of formation have been found for systems 
larger than a pair of qubits, except for some special symmetric states [15]. 

For a multipartite quantum system, the degree of entanglement will neither increase 
nor decrease under local unitary transformations on a quantum subsystem. Therefore the 
measure of entanglement must be an invariant of local unitary transformations. The entan- 
glements have been studied in the view of this kind of invariants and a generalized formula 
of concurrence for high dimensional bipartite and multipartite systems is derived from the 
relations among these invariants [16]. The generalized concurrence can be used to deduce 
necessary and sufficient separability conditions for some high dimensional mixed states [17]. 
However in general the generalized concurrence is not a suitable measure for iV-dimensional 
bipartite quantum pure states, except for N = 2. Therefore it does not help in calculating 
the entanglement of formation for bipartite mixed states. 

Nevertheless in [18] it has been shown that for some class of quantum states with N > 2, 
the corresponding entanglement of formation is a monotonically increasing function of a gen- 
eralized concurrence, and the entanglement of formation can be also calculated analytically. 
Let Ti be an N- dimensional complex Hilbert space with orthonormal basis e*, % — 1, ...,N. 
A general bipartite pure state on Ti <8> Ti is of the form, 

N 

\i> >= a V C i ® € V a ij e G (!) 

N 

with normalization ^ a ij°*ij — 1- The entanglement of formation E is defined as the 
entropy of either of the two sub-Hilbert spaces [8], 

E(m = -Tr ( Pl log 2 Pl ) = -Tr (p 2 log 2 p 2 ) , (2) 

where p\ (resp. p 2 ) is the partial trace of over the first (resp. second) Hilbert space 

of Ti <8> T~t. Let A denote the matrix with entries given by in (1). p 1 can be expressed as 
Pl = AAl 
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The quantum mixed states are described by density matrices p on H®H, with pure-state 
decompositions, i.e., all ensembles of states of the form (1) with probabilities Pi > 0, 
p = Y J \=iPi\i>i)('4>i\- l Y,i=iPi = 1 for some / G IN. The entanglement of formation for the 
mixed state p is defined as the average entanglement of the pure states of the decomposition, 
minimized over all decompositions of p, E(p) = min Yh=i PiE(\4>i))- 

For N = 2 equation (2) can be written as E(\ip))\ N=2 = h((l + y/l - C 2 )/2), where 
h(x) = — x\og 2 x — (1 — x) log 2 (l — x). C is called concurrence, C{\ip)) = 2|ana 2 2 — 012021 1 
[14]. E is a monotonically increasing function of C and therefore C can be also taken as 
a kind of measure of entanglement. Calculating E(p) is then reduced to the calculation of 
the corresponding minimum of C{p) = min J2if=i PiC (\4>i)) , which simplifies the problems, 
as C{\ipi)) has a much simpler expression than E{\ipi)). 

For iV > 3, there is no such concurrence C in general. The concurrences discussed in 
[16] can be only used to judge whether a pure state is separable (or maximally entangled) 
or not [17]. The entanglement of formation is no longer a monotonically increasing function 
of these concurrences. 

Nevertheless, for a special class of quantum states such that AA^ has only two non-zero 
eigenvalues, a kind of generalized concurrence has been found to simplify the calculation of 
the corresponding entanglement of formation [18]. Let Ai (resp. A 2 ) be the two non-zero 
eigenvalues of AA^ with degeneracy n (resp. m), n + m < N, and D the maximal non-zero 
diagonal determinant, D = A™ A™. In this case the entanglement of formation of is 
given by E(\ip)) = —n\i log 2 A x — mA 2 log 2 A 2 . It is straightforward to show that E(\ip)) is 
a monotonically increasing function of D and hence D is a kind of measure of entanglement 
in this case. In particular for the case n — m > 1, we have 



d is defined to be the generalized concurrence in this case. Instead of calculating E(p) 





where 




and 




(4) 
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directly, one may calculate the minimum decomposition of D(p) or d(p) to simplify the 
calculations. In [18] a class of pure states (1) with the matrix A given by 



/ b ai 6i \ 

—6 Ci di 

ai ci — e 

V 6i di e / 



(5) 



ai, bi, Ci, di, 6, e G C, is considered. The matrix AA^ has two eigenvalues with degeneracy 
two, i.e., n = m = 2 and \AA^\ = \b\C\ — a\d\ + be\ A . The generalized concurrence d is 
given by d = 4|fe 1 Ci — ai^i + 6e|. Let p be a 16 x 16 matrix with only non-zero entries 

Pl,W = P2.15 = -P3.14 = P4.10 = P5.12 = P6,H = P7.13 = ~P8,8 = ~P9,9 = Pl0,4 = Pll,6 = Pl2,5 = 

Pi3 ; 7 = — Pi4,3 = Pi5,2 = Pi6,i = 1- can be further written as 

d = m P r)\- (6) 

Let * denote the set of pure states (1) with A given as the form of (5). Consider all mixed 
states with density matrix p such that its decompositions are of the form 

M M 

p = 5>^>(^|, & = hk)etf. (7) 

i=i i=i 

All other kind of decompositions, say decomposition from |^), can be obtained from a 
unitary linear combination of [14, 18]. As linear combinations of do not change 
the form of the corresponding matrices (5), once p has a decomposition with all \ipi) e VP, 
all other decompositions, including the minimum decomposition of the entanglement of 
formation, also satisfy that \ip[) G Then the minimum decomposition of the generalized 
concurrence is [18] 

16 

d(p)=A 1 -J2K, (8) 

i=2 

where Aj, in decreasing order, are the square roots of the eigenvalues of the Hermitian matrix 



R = yj y/ppp*Py/p, or, alternatively, the square roots of the eigenvalues of the non-Hermitian 
matrix ppp*p. 

2 Entanglement of formation for a class of high dimen- 
sional quantum states 



An important fact in obtaining the formula (8) is that the generalized concurrence d is a 
quadratic form of the entries of the matrix A, so that d can be expressed in the form of 
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(6) in terms of a suitable matrix p. Generalizing to the iV-dimensional case we call an pure 
state (1) (^-computable if A satisfies the following relations: 

\AAl\ = ([A][A]*)W 

9 

\AA* - XId N \ = (A 2 - \\A\\X + [A] [A]*) N/2 , 

where [A] and \\A\\ are quadratic forms of a^-, Idx is the N x N identity matrix. We denote 
A the set of matrices satisfying (9), which implies that for A e A, AA^ has at most two 
different eigenvalues, each one has order N/2 and d is a quadratic form of the entries of the 
matrix A. 

In the following we give a kind of constructions of high dimensional rf-computable states. 
For all N 2 x N 2 density matrices with decompositions on these iV-dimensional (i-computable 
pure states, their entanglement of formations can be calculated with a similar formula to 
(8) (see (31)). 

We first present a kind of construction for a class of TV-dimensional, N = 2 k , 2 < k E IN, 
d-computable states. Set 

A 2 = 

where a,c,d G C. For any b\,c\ e C, a 4 x 4 matrix A 4 <E A can be constructed in the 
following way, 




A 4 



Bo A, 



— a 1 r l 



where 



B 2 — b\J<ii C*2 — C\Jii J 2 



( b\ a —c 

-61 c d 

—a — c —c\ 

c —d Ci 



(10) 



/ 1 \ 
-1 



V 



/ 



A 4 satisfies the relations in (9): 



A A A\ 



*\2 



[(hd + ad + c^ihd + ad + c 2 )*] 2 = ([MIM*) 
A 4 A\ - \Id 4 = (A 2 - + c x c\ + aa* + 2cc* + dd*)\ 
+ (&1C1 + ad + c 2 )(6iCi + ad + c 2 )*) 2 
= {X i-\\A 4 \\\+[A 4 ][A 4 }*y, 
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where 



[A 4 ] = (61C1 + ad + c 2 ), ||At|| = + cic* + aa* + 2cc* + dd* 
A 8 G ^4 can be obtained from At, 

A s = 

where 

-B 4 = &2^4, C4 = c 2 J 4 , ^4 = , h, c 2 G <C. 



(11) 




(12) 



/ 1 \ 
10 
0-100 
-10 



(13) 



For general construction of high dimensional matrices A 2 k+i G A, 2 < k G IV, we have 

/ B 2 k A 2 k \ ( bk-J 2 k A 2 k ^ 

k(k+i) 



A 



2 k+i 



\ (-I)" 



2 A^k C^k 



y (-1) 2 A* fc Cfcj^ y 



(14) 



J 2 k + 1 



V 



:-i) 





(fc + l)(fc+2) 



J 2 fc \ 



j^J; 



(15) 



where b k , c k G <C, -B 2 fc = b k J 2 k, C 2 k = c k J 2 k. We call J 2 k+i multipliers. Before proving that 
A 2 k+i e we fi rs t give the following lemma. 



Lemma 1. A 2 k+i and J 2 fe+i satisfy the following relations: 

J 2 k + 1 J 2 k+1 = J 2 k+1 J^k + l — Id 2 k+1, 

J 2 k+iJ 2 k+i = J 2 k+iJ 2 k+i = (—1) 2 Id 2 k+1, 



"2 k + 1 — <- , 2 fe + 1 ' 



_ (fc+i)(fc+2) T 

J 2 fe + 1 — — 1J 2 J 2 k + 1, 



.hih+R 

2 /In 



(16) 



(17) 



A 2 k+i — {— 1) 2 -<4-2*+ 1 j A 2k+1 = (—1) ■= j-L 2 k+i- 
Proof. One easily checks that relations in (16) hold for k = 1. Suppose (16) hold for 
general k. We have 

I (-i) ih±1 ¥ m jt k \f J 2k \ 

J 2 k 



J 2 k + lJ 2 k+l — 



\ 



) 



(fc+l)(fc+2) 

-1) 2 J* fc 



V 



/ (-l)(*+l)(*+2)j 2fc ^ \ 



V 



Id 2 k + \ 



J 



and 



Tt Jt 

d 2 k+i ^ 2 fc+l 



( (-1)^^ \ 



J 2 k 

/ / ,s (fc + l)(fc+2) j Jt 

(-1) 2 Jrffc 



( (-1)^^^ \ 



J 2 fc 





(fc + l)(fc+2) 



V 



o (-1) — 2 — J l 2k J 2 k 



Therefore the relations for J* fc+1 J 2 fc+i and J 2k+1 J 2k+1 are valid also for /c + 1. The cases for 
J 2 fc+i J^+i an d ^2 fc + 1 ^2 fc + 1 can be similarly treated. 

(fc+1) (fc + 2) 

The formula J 2fc+ i = (— 1) 2 J 2 i=+± in (17) is easily deduced from (16) and the fact 

Tt _ jt 

The last two formulae in (17) are easily verified for k — 1. If it holds for general k, we 
have then, 



/ 



fc(fc+i) 



B\ k (-l)^s- A 2 h 



A 1 



C 2 fc 



I k(k+l) fc(fc+l) . \ 

/ ("I) 2 ^ 2fc (-1) 2 A 2 fe \ , fc(fc + l) 

= (-1) 2 ^4 2 fc+i, 

i.e., it holds also for + The last equality in (17) is obtained from the conjugate of the 
formula above. ■ 



Lemma 2. The following relations can be verified straightforwardly from Lemma 1, 



. fc(fc+i) fc(fc+i) 

1) 2 B 2 fc, C* fc = (-1) 2 C 2 fe, 



2 fc — V — ± J * J -'2 K i < - / 2 fc 

D-l — 1 R< 1 pt ^-1 _ 1 /nrt _ _J_r»t 

-°2fe+i ~~ l2- D 2 fe + 1 — 717 2 k +ii <^ 2 fc+i — — 2 fc+i — -0 2fc+ i. 

b k b k b k 



CkC k 



B 2k+1 B 2 k+i — B 2 k+iB 2k+1 — b\ld 2 k+i, C 2 k+iC 2 k+i — C 2 k+iC 2k+ i — c\Id 2 k+i, 
B 2k+1 B 2 k+i = B 2 k+iB 2k+1 = b k b^I d 2 k+i , C 2k+1 C 2k +i = C 2 k+iC 2k+1 = c k c* k Id 2k +i. 

For any A 2 k+i E A, k > 2, we define 

||A 2 fc+i|| =: b k b k + c k c k + | \A 2 k ]|, 
[A 2fc+1 ] =: (-l)^+D/2 6fcCfe _[A 2fe ]. 

Lemma 3. For any k > 2, we have, 

(A^+i ,/ 2 fc+i ) (,/ 2 fc+i A 2 fc+i)' = (A 2 fc+i J 2 fe+i)'(^ 2 fc+iA 2 /c+i) 

k ( fc~(-i) 

= ((-l)^ - &fcc fe - [A 2 fc])i"d 2 fc+i = [A 2 fc+i]/d 2 fc+i, 
(A* fe+ i J 2 fc+i)(J 2 fc+iA* fc+ i)' = (A* k+1 J 2 fe+i)*( J 2 fc+i A* fe+ i) = [A 2 fc+i]*7d 2 fc+i. 



(18) 

(19) 
(20) 



(21) 



(22) 



Proof. One can verify that Lemma 3 holds for k = 2. Suppose it is valid for k, we have 

( A 2 fc+l J 2 k+l) (<^2 fc + 1 A 2 k+l Y 



it 

J 2 k 



(fc+l)(fc+2) fc(fc+l) / 1 - (fc + l)(fc+2) , (fc+l)(fc+2) 



en e i2 
e2i e 2 2 

where 



(fc+l)(fc+2) + fc(fc + l) fc(fc+l) 

en = (-1) 2 A 2 fc J t 2k A 2 kJ t 2k + (-1) 2 b k c k Id 2 k 

(fc+l)(fc+2) + fc(fc-l) , , , fctfc+1) 

= (-1) 2 (A 2fc ^ fc )(' / 2*^2*) + (-l)^^b k c k Id 2k 

= ((-lJ^^ftfcCfc - [A 2 fc])/tZ 2 fc, 
ei2 = b k A 2 k J\ k + (-1) 2 b k A\ k J 2 k 

U A Tt /1 , / . (fc+l)(fc+2) + fc(fc-l) , 

= 6 fc A 2fc J^ fc (l + (-1) 2 )=o, 
e 2 i = (-l) (fc+1 2 fc+2) c fc A 2 fc J* fc + (-l)^^ 11 ^^ J 2 fc = 0, 

fc(fc+l) , r , , (fc+l)(fc+2) + fc(fc+l) f T A t T 

e 2 2 = (-1) 2 b k c k Id 2 k + (-1) 2 A^ fc J 2 fcA^J 2 fc 

fc(fc+l) (fc+l)(fc+2) + fc(fc-l) 

= (-1) 2 b k c k Id 2 k + (-1) 2 (A 2 /= J^XJ^A^J* 

= ((-l)^~hc k - [A 2 k\)Id 2 k. 



Hence 

fc(fc+i) 

( y A 2 k+iJ 2 k+i)( y J 2 k+iA 2 k+i') = ((—1) 2 b k c k — [A 2 k])Id 2 k+i = [A 2 k+i])Id 2 k+i. 

Similar calculations apply to (A 2 k+i J 2 k+if{ J 2 k+iA 2 k+i). Therefore the Lemma holds for k + 1. 
The last equation can be deduced from the first one. 

Theorem 2. A 2 k satisfies the following relation: 

\A 2k+i A\ k+1 \ = ([A 2k+1 ][A 2 k +1 ]*f = [((-l)^b k c k - [A.kM-l^btd - [A 2 k]*)f. 

(23) 
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Proof. By using Lemma 1-3, we have 

^ B 2 k A 2 k ^ 

fctfc+1) , 

V (-1)^-^*2* 



|A 2 fc+i| = 



Therefore 



C\ k f 







J 2 h > 

Id 2 k 

fc(fc+i) 



5 



2 fc 



A 2k \ 



2 4* C\ k ) 



B 2 k - {-\)^ A 2 k{C\ k y l A\ k o \ 



A 2 k (C 2k/ 
fc(fc+i) 



fc(fc+i) ^ 

\b k c k Id 2 k - (-l)^ 1 ^ (A 2 kJ 2 k)(J 2 kA 2k y\ 

b k c k Id 2 k - [A 2 k]Id 2 k \ = ((-l) k ^ b k c k - [A 2 fc]) 2fc . 



-1) 2 



|A 2 fc+i A!, fc+ i | = ([A 2 fe+i][A 2 fc+i]*) 2 . 



Lemma 4. A 2 fc+i and J 2 fc+i satisfy the following relations: 

(A 2 k+i J 2 k+i)( J 2 k+iA 2 k+i)^ + ( J 2 k+iA* 2k+1 ){ J 2 k+iA 2 k+iY 

= A 2 k+iA 2k+1 + J 2 k+i A 2k+1 A^k+i J 2k +i = 1 1 A 2 k+i | \Id 2 k+i , 

(A 2 k+iJ 2 k+iy(A^ k+1 J 2 k+i) + ( J 2 k+i A 2 k+iy( J 2 k+i A 2 k+i) 

= A 2k+1 A 2 k+i + J 2k+1 A l 2k+1 A* 2k+1 J 2 k+i = \ \A 2 k+i\\Id 2 k+i. 

Proof. It can be verified that the first formula holds for k — 2, if it holds for k, we have 
(A 2 k+i J 2 k+i)(A 2 k+i J 2 fc+i)^ + ( J 2 fc+iA2 fe+ i)( J 2 fc+i A 2 fc+i)' 



/ (fc+l)(fc+2) 
(-1) 2 



^4 2 fe J^fc 



fc(fc+i) 



fc(fc+i) 

"1) 2 ^4 2 fc<^2 fe 



2 fe 



+ 



W-i) 

/ll 7l2 
/21 ^22 



1) 2 

(fc + l)(fc+2) 



^2 fc ^ 2 fe ^2 fe C 2 fc 
J 2 kB* k ( — 1) 



(fc + 1) 2 (fc+2) T* 4* / 
2 J 2fc j A l 2 fc / 



V 

V 



2 fc 2 fc 



l)^^A 2fc J* fc (-1) 2 B\ k J 2 k 

(fc+l)(fc+2) 



"I) 



fe(fe+l) 
- , ; 7* A* 



'2 fc 2 fe 



(fc + l)(fc+2) 



C 2 k J 2 k 
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where, by using Lemma 1 and 2, 

/n = / 22 = A 2k A\ k + J 2k A\ k A 2k Jl k +BB^ + J 2 k&CJ* k 

j. fc(fc+i) ^(fc+i) . . 

= A 2k A\ k + J 2fc (-1)— 2— A£ fc (-1)— 3— 4** J\ k + (& fc &£ + c k c* k )Id 2 k 

= (\\A 2 k\ \ + b k b* k + c k c* k )Id 2k = \\A 2k+ i\\Id 2 k, 

X A ^y* , f s fc(fc+l) fc(fc + l) + (fc+l)(fc+2) | (fc+l)(fc+2) 

/12 = A 2 fcC* fc + (-l) 2 £ 2fc A; fc + (-l) 2 6 fe J 2fc ^ fc + (-l) 2 c* k A\ k J 2 k 

= ( -i)^(^^ + (-i) fc(fc - 1)+ r Kfc+2) ^^) 

+A 2fc q fc + (-i) fc(fc - 1)+( r )(fc+2) A 2fc c 2 ^o, 

/ 21 = c 2 »4» + (-i)^^^ + (-i) fc(fc+1)+( 2 +1)(fc+2) ^ 2fc ^ fc + (-1)^.41 



2 fc 

- fc(fc+l) , t t , fc(fc-l) + (fc+l)(fc+2) t t . 

= (-1) 2 (6 fe A 2fc J* fe + (-1) 2 6 fc A 2fc J* fc ) 

+c^J 2fc A 2fc + (—1) 2 c£J 2fe A 2fc = 0. 
Hence the first formula holds also for k + 1. The second formula can be verified similarly. ■ 

Lemma 5. Matrices B 2 k, A 2 k and C 2 k satisfy the following relations: 

k(k-{-l) k(k-\-l) 

({-l)^B 2k A* 2k + A 2k C* k ){{-\)^- A\ k B 2k + C; k A 2k y = F(A 2k+ i)Id 2k , (24) 

where 

k(k-\-l) 

F(A 2k+1 ) = cf[A 2k ] + b 2 k [A 2k ]* + (-l)—b k ci\\A 2k \\. (25) 



Proof. By using Lemma 3 and 4, we have 



fc(fc+i) . fc(fc+i) 



((-i)^s 2fc ^ t + A 2fc c* fc )((-i)^^*5 2fc + c* k A 2k y 



— b k (J 2k A 2k )(A 2k J 2k y + c k f(A 2 kJ 2 k)(J 2 kA 2 k) 



+ (_i)^6 fcC *[(^ fc j 2fc )(A; fc j 2fc ) 4 + (j 2fc ^ fc )(W 2 *)*] 



= (c£ 2 L4 2 *] + bt[A 2 k]* + (-1)^6^11^11)/^ = F(A 2fc+1 )7d 2fe . 
Lemma 6. A 2 k and J 2 fc satisfy the following relation: 

\\A 2 k || J 2 kA* 2k A\k J 2 k = [A 2 k] [A 2 k]*Id 2 k + J 2k A\ k A\ k A* 2k A\ k J 2k . (26) 
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Proof. From (24) we have the following relation: 

F( A 2 k+i ) J 2 k A 2k A 2k J 2k 

= ((-l)^B 2 kA* 2k + A 2k C* k )(-l) ! ^b k (A* k J 2 ,) t J 2k A; k A\ k J t 2k 

+(( — 1) 2 B 2 kA* k + A 2 kC 2k )c%{J 2 kA 2 k) t J 2 kA* 2k A t 2k J 2k 
= (-l) kJ ^ bk ((-l) kJ! ¥ 1 B 2k A* k + A 2fc C 2 \)[(^ fc J 2fc )*(J 2 ^; fc )]^^ 

-0 2 -^2 fc ^2 fc (^2 fc ^2 fe )* t -^2' c A 2k A 2k J 2k ^ A 2 kC 2 k(J 2 kA 2 kY J 2 kA 2k A 2k J 2k ~\ 
= (-l)^b k ((-l) k ^b k J 2k A* k [A 2k ]Al k Jl k +c* k A 2k J 2k [A 2k }*A t 2k J t 2k ) 

~ >t ~ c k[(~ -0 2 b k J 2 kA 2k A 2k A 2k A 2k J 2k +c k [A 2k ]J 2 kA 2k A 2k J 2k \ 
= b 2 k [A 2k }* J^A^A^J^ + (-l) ! ^b k d' k [A 2 k][A 2 k]*Id 2 k) 

+(—1) 2 c k b k J 2k A 2k A 2k A 2k A 2k J 2k +c 2 [A 2k ]J 2k A 2 kA 2k J 2k ] 

= (b 2 k + cD[A 2 k]J2kA 2 kA t 2k Jl k + (-l)* I T tli 6fcC fc ([A 2fc ] 2 /d 2 * 

fc(fc+i) 

+(—1) 2 c k b k J 2 kA 2k A 2k A 2 kA 2k J 2k ). 
Using (25) we have 

11^2* II J2 k ^ 2 kA 2k J 2 k = [A 2 k][A 2 k]* I d 2 k + J 2 kA 2k A 2k A 2k A 2k J 2k . 



Theorem 3. The eigenvalue polynom of A 2 k+iA' 2k+1 satisfies the following relations 

\A 2k+1 A\ k+1 - XId 2k+1 \ = (A 2 - P 2fc+ i||A + [A 2k+1 }[A 2 k +1 ]*) 2k , 
\A\ k+1 A 2 k +1 - A/rf 2fc+ i| = (A 2 - p 2fc+ i||A + [A 2 k +1 ][A 2k+1 ]*) 2k . 



Proof. Let 



, k(k+l) -i 

Afe = -[(cfcCfc - A)/rf 2fc + ^ fc ^ fc ][(-l)— a - B 2k A* 2k + A 2fc C 2 * fc ] _1 . 
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A 2 k+lA 2k+1 — \Id 2 k + l 




/ (-l) hih T 11 B 2k A* 2k + A 2k C* k (b k b* k - \)Id 2k + A 2 uA\ k \ 
(cfcCfe ~~ X)Id 2 k + A 2k A 2k (— 1) 2 A 2k B 2k + C 2k A 2k 

/ (-\) hih f R B 2 kAl k + A 2 kC* 2k {b k bl-\)Id 2 k+ A 2 kA\ k ^ 
( c fc c fc ~~ X)Id 2 k + A 2k A 2k (— 1) 2 A 2k B 2k + C 2k A 2k 

(b k b* k - X)Id 2 k + A 2 kA\ k 



-Vf^ 1 B 2 kA* 2k + A 2 kC, 



\ 







= \I + U\, 
where 

/ = ((-i)^s 2fc ^ fc + A 2fc C* fc )((-l)^S* fc A 2 * + A; fc c 2fc )' 

fc(fc + l) k(k-\-l) x 

= (-l)^6 fc c fc [A 2fc ]*/rf 2fc + (-l)^b* k c* k [A 2 k]Id 2k + b k b* k J 2 kA* 2k A\ k J\ k + c fc ^A 2fc ^5 
and, by using Lemma 5, 

// = -{{-xf^B^AI^ + A 2fc C* fc )A fc [(& fc &£ - \)Id 2k + A 2 *A 
= [(c k c* k -A)((-l) 



"2* 



l 2 fc] 



2 i^fcA^fc + A^C^fc 



^2 fc C 2 fc) "I - (( 1) 2 -B2 fc ^. 2 fc "I" A 2k C 2k )A 2k A 2k ] 
[(b k b* k - A)((-l) M ^ ii J B 2fc ^ + a^c*,)" 1 + ((-l)^B 2k A* 2k + A 2fc c* fc ) _1 4^ 

- A)(c„4 - A)/rf 2fc + (b k b* k - X)((~l) kJ ^B 2 kA; k 
-\-A 2k C 2k )A 2k A 2k (( 1) 2 .B^AJjfc + A 2 kC 2k ) 1 + (c k c k A)A 2 fcA2 fc 
~K( — 1) 2 B 2 kA 2k + A 2 kC 2k )A 2k A 2k ((— 1) 2 B 2 kA* 2k + A 2k C 2k ) 1 A 2k A 2k 
(b k b* k - \)(c k c* k - X)Id 2 k + (c fc c* - A)A 2fc A 2fc + rtX^T) 111 + F(1^7) J^2*4*> 
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where 

/// = ((-i)^s 2fc ^ fc + A 2k c* 2k )A t 2k A* k ((-i) m ^ 1 A* 2k B 2k + c; k A 2k y 

fc(fc+l) ^ ^ . . * // \ fc ( fc + 1 ) * * \t 

= ((—1) 2 B 2 kA 2k + A 2 kC 2k )A 2k J 2k J 2 kA 2k {{— 1) 2 A 2k B 2 k + C 2k A 2 k) 

k(k-\-l) 

= [(-l)^b k (J 2 kA* 2k )(J 2 kA 2 kY + c* k (A 2 kJ 2 k)(J 2k A 2 ky] 

k(k-\-l) 

[{-l)^b k (J 2 kA* k )(A* k J 2 kY + c* k {J 2 kA* k )(J 2 kA 2 ky] 

k(k-\-l) 

= [(-l)^b k (J 2 kA* k )(J 2 kA 2 ky + c%[A 2k \Id 2 k\- 
[(-l)^6 fc [^ 2fc ]*/d 2fc + ^(J^X^Aj*)*] 

= (6^[A 2 fc]* + C^ 2 [A2fc]) t /2fcA2fc A^kJ^k + ( — 1) 2 bkC^^kA^kA^kA^kA^kJ^k 
fcffc+i) 

From Lemma 6, we get 

/// = (bl[A 2 k]* + cl 2 [A 2 k])J 2 kA; k Al k Jt k + (-l) k ^ 
= F(A 2k +i)J 2k A 2k A 2k J 2k . 

From Lemma 3 we also have 
IIIA 2 kA 2k = IIIA 2 kJ 2 kJ 2 kA 2k 

— F(A 2 k+i)J 2 kA 2k (J 2 kA 2 k) t (A 2 kJ 2 k)(A* k J 2 k) t = F(A 2k +i)[A 2 k][A 2 k]*Id 2k . 

Therefore, 



A 2k +iA 2k+1 — XId 2 k+i 



\I + H\ 



k(k-\-l) 

= I - A 2 /rf 2fc + \{b k b* k + c k c* k + \\A 2k \\)Id 2 k - (b k b* k c k c* k - {-l)—b* k c%[A 2 k\ 

k(k-\- 1) 

- (-l) — b k c k [A 2k ]* + [A 2 k][A 2 k]*)Id 2k \ 
= (X 2 -\\A 2k+1 \\X+[A 2 k +1 ][A 2 k +1 Y) 2k , 
where the first formula in Lemma 4 is used. The second formula in Theorem 3 is obtained 
from the fact that A 2k +iA 2k+1 and A 2k+1 A 2k +i have the same eigenvalue set. ■ 

From Theorem 2 and 3 the states given by (14) are ^-computable. In terms of (4) the 
generalized concurrence for these states is given by 

d 2k+ i = 2 k+1 \[A 2 k+i} \ = 2 k+1 \b k c k + 6fe-ic fe _i + ... + bid + ad + c 2 \. 

Let p 2 k+i be a symmetric anti-diagonal 2 2k+2 x 2 2k+2 matrix with all the anti-diagonal 
elements 1 except for those at rows 2 k+1 - 1 + s(2 k+2 - 2), 2 k+1 + s(2 k+2 - 2), 2 k+2 - 1 + 
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s (2 fc + 2 _ 2), 2 k+2 + s{2 k+2 - 2), s = 0, 2 k+1 - 1, which are -1. d 2k +i can then be written 

as 

d 2 k+l = \{lp2 k + 1 \P2 k + 1 i J 2 k + 1 )\i (28) 

where 

2 fc+i 

l^+i) = XI (^2 fc+1 )«i e * ® e i- ( 29 ) 

Let $ denote the set of pure states with the form (29). For mixed states with density- 
matrices such that their decompositions are of the form 

M M 

p 2 2k +2 = J2pM(A\, = i, \A) e (30) 

i=l i=l 

their entanglement of formations, by using a similar calculation in obtaining formula (8) 
[18], are then given by E(d 2 k+i (p2 2fc + 2 )), where 

2 2k+2 

d 2 k+i (p 2 2k+2) = Qi — X (31) 

and fij, in decreasing order, are the the square roots of the eigenvalues of the matrix 
p 2 2k+2p 2 k+i p* 2k+2 p 2 k+i . Here again due to the form of the so constructed matrix A 2 k+i in 
(14), once p has a decomposition with all G $, all other decompositions of |^) also 
satisfy |^) G <&. Therefore from high dimensional ci-computable states A 2 k+i, 2 < k < N, 
the entanglement of formation for a class of density matrices whose decompositions lie in 
these ci-computable quantum states can be obtained analytically. 

3 Remarks and conclusions 

Besides the d-computable states constructed above, from (10) we can also construct another 
class of high dimensional ci-computable states given by 2 k+1 x 2 k+1 matrices A 2 k+i , 2 < k G 
N, 

( B k A k \ ( hkhk A2k 



A 2 k + 1 = 

where bk, c k G C, h — Ja, 



v ~ A k C k J y _ A t^ Ck j 2k 
I l 2 k\ 



(32) 



I 2 k+1 



-I 2k 

V / 



(33) 
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for k + 2 mode 4 = 0, 
for k + 2 mode 4 = 1, 

I 2 k+l = 

for /c + 2 mode 4 = 2, and 

/ 











I 2 k-2 

\ ~I 2 k-2 

for k + 2 morfe 4 = 3. 



J 2 fc+i 



I 2 k+l 


-(,° 


I 2 k \ 
o / 













I 2 k-1 








—I 2 k- 





n 

u 


T t i 
i 2 fc-l 


n 

u 


n 

u 


-l 2 k- 


n 

-1 u 


n 

u 


n 

u 


(J 


(J 


(J 


(J 


(J 


(J 


(J 


u 


u 


u 


u 


— i 2 fc-2 








T 

l 2 k-2 








— I 2 k-2 








J 2 fc-2 




































1 \ 





— I 2 k-2 










I 2 k~2 ^ 








) 



One can prove that the matrices in (32) also give rise to ^-computable states: 

k k 



\*l2* 



|A 2 fc+i^ fc+1 | = [(c 2 + ad - ^bid)(c 2 + ad - ^hd)* 

i=i i=i 

k k 

\A 2k +iA\ k+1 - \Id 2 k+i | = [A 2 - (aa* + 2cc* + dd* + ^ hb* + ^ qc* ) A 



i=l i=l 
k 



(34) 



(35) 



(36) 



+ ( c 2 + arf _ J2b iCi )(c 2 + ad -J2kc t yf- 

i=l i=l 

The entanglement of formation for a density matrix with decompositions in these states is 
also given by a formula of the form (31). 

In addition, the results obtained above may be used to solve linear equation systems, 
e.g., in the analysis of data bank, described by Ax = y, where A is a 2 k x 2 h matrix, k G IN, 
x and y are 2 fc -dimensional column vectors. When the dimension 2 k is large, the standard 
methods such as Gauss elimination to solve Ax. = y could be not efficient. From our Lemma 
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3, if the matrix A is of one of the following forms: A 2 k, B 2 kA 2 k, A\ k or A^B^, the solution 
x can be obtained easily by applying the matrix multiplicators. For example, A 2 fcX = y is 
solved by 

x = — — - (A 2 k J 2 k) 1 J 2*y- 

[A 2 k\ 

The solution to B 2 kA 2 kK = y is given by 

X = U \A 1 (^2*^2*)* J 2 fc Y- 

bk[A 2 k\ 

We have presented a kind of construction for a class of special matrices with at most 
two different eigenvalues. This class of matrices defines a special kind of ^-computable 
states. The entanglement of formation for these ^-computable states is a monotonically 
increasing function of a the generalized concurrence. From this generalized concurrence the 
entanglement of formation for a large class of density matrices whose decompositions lie in 
these ci-computable quantum states is obtained analytically. Besides the relations to the 
quantum entanglement, the construction of rf-computable states has its own mathematical 
interests. 
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